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Abstract 

Motivated by recent developments on cold atom traps and high density QCD we consider 
fermionic systems composed of two particle species with different densities. We argue that a 
mixed phase composed of normal and superfluid components is the energetically favored ground 
state. We suggest how this phase separation can be used as a probe of fermion superfluidity in 
atomic traps. 
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The recent interest in two physical systems have revived the study of asymmetric 
fermionic systems, that is, systems with unequal number density (or chemical potential) 
for the different species. In high density strongly interacting systems, as it may be found at 
the core of "neutron" stars, the different quark flavors have different chemical potentials on 
account of their different masses and charges Q3H 3 Atomic traps can also provide an 
example of similar systems if a bias can be introduced when the trap is filled in order to have 
a larger number of one of the atom species (or hyperfine state of the same atom) (Hf . In the 
symmetrical situation the low energy properties of these systems are dominated by Cooper 
pair formation. Since it involves the attraction between two fermions with equal and oppo- 
site momenta at the their Fermi surface one could imagine that as the Fermi surfaces move 
apart with increasing asymmetry, the pairing would become weaker and the gap smaller, 
until superfluidity disappears. What actually happens, the formation of an inhomogenious 
mixed phase state, is however more interesting and results from the competition between 
states with different particle distribution, both in momentum and real space. 

Our discussion is valid in a wider class of models but, for definiteness, let us consider a 
non-relativistic dilute gas made out of two particle species A and B with chemical potentials 
Ha, Hb and masses Ma, Mb respectively. At low densities the details of the potential are 
not probed and their interaction is well described by the pairing Hamiltonian 
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where ip A ,ipA are creation and annihilation operators for the A particles and is their 
dispersion relation, that we take to be = k 2 /2M A — \ia (and similarly for B). In the 
mean field approximation, adequate for the low densities considered here, the Hamiltonian 
can be approximated by 
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where A = -g J d 3 kj '(2n) 3 {4j B (-k)^ A {.k)) 
ef)/2 and the fields ip a ,ip^ are defined by 
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It is straightforward to minimize the diagonalized Hamiltonian shown in Eq. (J2J) . One simply 
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fills the modes with negative Eu and leave the remaining modes empty. More precisely, 
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FIG. 1: Dispersion relation for the quasi-particles a and (3 showing a region where E% is negative 
for M B = 7Ma, Pb = 1-^pa, and A = 0.02/i^ < (p 2 B - p a )/(A^/MaM b ). Solid curve correspond 
to E^ (top) and E^ (bottom). 



the ground state satisfies 

V>«,„(A;)|*) = if Ef>0, 

In terms of the original particles A and B and the vacuum state |0), the state above 
corresponds to having a BCS-like state [u^ + Vkip A {k)il> B (—k)}\0) in the modes k where 
E^ > 0, but a state filled with particle B (A) in the modes where E% < (E% < 0). The 
thermodynamic potential of this state is 

'6(-E£)E£ + 0(-EZ)EZ + 4 - El , 
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where a is the scattering length between particle A and B. It is related to the coupling 
constant (in dimensional regularization) by 1/g = M/(2ira) with M = MaMb j (Ma + Mb), 
the reduced mass. 

Let us now consider the case M# > , ps > Pa, where Pi is the Fermi momentum 
defined by pi = ^j2M^Xi. Fig. [T] shows that, for some values of A, E® may be negative for 
momenta k± < k < k 2 where 
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k\ 2 = V - A ^ ± g y/iPk - - 16M A M B A2 , (7) 

while .E^ is always positive. 

We now discuss separately the cases where either the chemical potentials or the densities 
of each species are kept fixed. 
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1. Fixed chemical potentials 

In Fig. |21 we show the thermodynamic potential as a function of A for different values 
of pa and ps, keeping the combination p\jM = p\/Ma + p b /Mb fixed, computed from a 
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FIG. 2: Thermodynamic potential for different values of pb and pa (constant po). The top curve 
corresponds to pa = Pb and the lower curves correspond to increasing values of \p B — p A \. 

numerical evaluation of Eq. JBJ). An analytical expression valid for A <C ^a,b is available 
but it is not very enlightening. 

For large enough A, k\ 2 are not real, E^ is always positive and the thermodynamic 
potential is unchanged from the ps = Pa case. However, for A < (p B —p 2 A ) / (4^MaMb), k\ 2 
are positive and the thermodynamic potential can be lowered by filling the states between k\ 
and k 2 with /3-type quasi-particles. The A = state, in particular, has its thermodynamic 
potential lowered with increasing p\ — p\ and at some point becomes smaller than the 
previous minimum with A = A corresponding to the BCS phase. The result is that, for 
fixed po and increasing p%—p\, there is a first-order phase transition between the superfluid 
and the normal state, as it has been noticed in different physics contexts (see, for instance 
Refs. Mil). These results can be understood in very simple physical terms. Suppose 
we start in the BCS phase with ha = A*b and increase ps- By absorbing one B particle and 
eliminating one A the system by one hand reduces its thermodynamic potential Q due to 
the —fiAfiA — bi^b term but by the other hand it increases Q by destroying two pairs. This 
is energetically favorable if, and only if, the difference in chemical potentials is large enough 
(or the gap small enough). Until that point the BCS state with equal number of particles 
remains the ground state, unchanged despite the variations in chemical potential. 

In addition to the stable (or metastable) normal and BCS phases there is, for some values 
of the chemical potentials fiA an d /i^, an unstable phase (referred to here as "Sarma state") 
corresponding to a maximum of Q as a function of A situated between the BCS minimum 
A and the normal phase at A = (first pointed out in Ref. The combination 

of parameters necessary for the existence of this phase can be found considering the gap 
equation: 




where the BCS gap is given by 



in the mean field approximation. For small values of the gaps A , A s <C fi A , the integrals 
can be approximated and it is found that 
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The Sarma state gap is in the range < A s < Ao, where the upper bound is set by the 
condition for the existence of real values of k\ 2 . We then have 

0<Ja ( - A ) < Ao =► Ao < < 2A . (11) 
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The condition in Eq. (fTT|) . for fixed p A and p B , can thus be seen as determining a window 
for the values of a supporting the Sarma state. It does not exist, even as an unstable state, 
for too small or too large interactions 0]. Also notice that when p 2 B — p\ reaches its largest 
value allowing for a Sarma state (p 2 B — p\ = 4a/M^Mb A ) , ki equals k 2 and the difference 
in particle densities n B — n A approaches zero, corresponding to a BCS state. In the other 
limit, p 2 B — p\ = 2a/ M a M b A , the gap A s vanishes and the Sarma state reduces to the 
normal phase. 



2. Fixed particle number 

The previous discussion regarding the stability of different phases was made under the 
assumption that the particle numbers n A and n B are allowed to change. We discuss now 
the situation where they are fixed. In the BCS phase the particle number densities n A and 
n B are the same, as can be readily seen by taking derivatives of Eq. (jUJ) in relation to fi A 
and n B , while the particle numbers in the normal phase can be different. The Sarma phase 
can also accommodate n A ^ n B so the question arises: what is the ground state of the 
system when the particle densities are fixed and different from each other? This question 
has been revived in Refs. [TH [3 EH where it was argued that the Sarma phase (named 
there "internal gap" ' state) could be the ground state in the case of fixed particle numbers 
n A , n B . This state would have fascinating properties, being at the same time a Fermi liquid 
(with two Fermi surfaces corresponding to k\ and k 2 ) and a superfluid. 

The question can be answered by finding the state with the smallest energy (not the 
thermodynamic potential Q). We compare here the the normal, Sarma and a mixed inho- 
mogeneous phase composed of bubbles of normal phase in a sea of the BCS phase. The 
energy for the normal and the BCS phase are, for small values of the gap, given by 
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A similar expression can be derived for the energy of the Sarma phase but its form is not 
very enlightening and it will not be needed below. 

The mixed phase is an inhomogeneous phase where a fraction x of the space is in the 
normal phase with A and B particle densities equal to n A and n B , while the remaining 
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FIG. 3: Figure on the left: Energy of the Sarma (dashed curve) and mixed phase (solid curve) 
relative to the normal phase as a function of 5n/n where n = (ha + ^b)/2 with Mb = 7 Ma- 
The values of ha and ns were obtained from the Sarma phase for a fixed a, and \a\pA = 0.59, 
0.63 < \a\pB < 0.65. Figure on the right: The fraction of normal state x as a function of bnjn 
with the same masses and scattering length. 



1 — x fraction is in the BCS phase with a common density for both species equal to n. The 
densities in each component are adjusted in such a way that the overall average densities 
have given prescribed values riA,nB, that is ua = xtia + (1 — x)n and similarly for the 
particles B. The most favored mixed state for given ua^b is the one with the smallest 
energy: 
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We have disregarded the interface energy between the two components, as those are small 
for large enough systems. Also, we assumed that the local asymmetries in the densities do 
not cost additional energy due to long range forces, as it would be the case if A and B had 
different charges 0]. We have to minimize Emix with respect to normal phase fraction x 
and the density in the BCS fraction n. 

There are two limiting cases where the comparison between the mixed and Sarma phases 
can be done analytically, corresponding to parameters where the inequalities in Eq. (|11|) are 
saturated. If p 2 B — p\ = 2\J 'MaMbA (A s =0), the Sarma state reduces to the normal state. 
In this case its energy is given by 
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where d~n = ris — ua is assumed to be small, 5n -C ua- 
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An upper bound on Emix — Es can be obtained by setting the density of the BCS 
component of the mixed phase n = ua and minimizing in relation of x. We have: 
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whose minimization yields 



(67T% A ) 4 /3 ^ n 2 



(16) 



Numerical calculations show that the upper bound above is close to the actual minimum. 
Notice that 
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thus, using A (p A ) = A (p ) + Q(Ao), we have 
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for Mb > M4 as it should be the case. Eq. (JT3J) shows that the mixed phase is energetically 
favored compared to the Sarma phase in one extreme of the window in Eq. (|11|). 

The other simple limit to analyze corresponds to P% ~ Pa = ^VMaMbA , in which case 
A s = A , ki = k 2 ., ua = n>B and the Sarma phase reduces to the BCS phase. In this case 
the mixed phase reduces to the BCS phase too and the energies of both the Sarma and 
mixed phase are equal to each other. For intermediate values of p 2 B —p\ ( still satisfying the 
constraint in Eq. ([lip), the difference Emix — Es interpolates between these two extremes, 
as Fig. (J3J) exemplifies. We find that for all reasonable values of the parameters (that is, 
where the mean field analysis should apply), and for fixed particle numbers via and ris, the 
mixed phase has a smaller energy than the Sarma phase. 



3. Fixed total density 

Another interesting situation arises when the total number of particles is fixed, but con- 
versions between particles of types A and B are allowed. This is relevant for the physics of 
high density quark matter where weak interactions can change the flavor of the quarks. In 
this situation n = ua + ub and 5fi = {ii B — A*a)/2 are fixed and the thermodynamic function 
that should be minimized is E — 5p,(n B — tia)- The non-relativistic formulation presented 
here is more appropriate for cold atoms. In cold atom traps, only species with nearly equal 
masses can convert into one another. Further, we consider Sfi = which is relevent for cold 
atoms for convenience. It is straightforward to see, then, that the condition for the existence 
of Sarma state Eq. (jllj) is not satisfied. The same also holds when 8/1 = and Mb 3> Ma- 
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The composition in the mixed phase is easy to calculate when 5/jl = 0. The most favored 
composition minimizes 



E(n = n A + n B ) = Min™ { (1 - x 
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When M^ = Mb, BCS (with x = 0) is the favored state and when Mb 3> M4 normal phase 
(sc = 1) with only particles of species B is favored. 

We have considered Fermi gases made up of two species, when an asymmetry on their 
densities or chemical potentials tries to push their Fermi surfaces apart, making pairing 
more difficult. We find that with either both chemical potentials or both densities fixed, the 
most likely ground state is a mixed phase composed of bubbles of an asymmetric normal 
state immersed in a sea of the symmetric BCS phase. It is worth to mention that the LOFF 
state (in which the condensate varies in space) has lower free energy than the normal 
and the BCS states. However, the LOFF state can exist only in a very narrow window of 
asymmetry for the chemical potentials and we ignored this posibility in our discussion. 

The space segregation of the excess particles in the mixed phase suggests a possible way 
of detecting superfluidity in atomic traps, specially where large gaps are expected as in 
the case of "resonance superfluidity" [17(. If an optical trap can be filled with an excess 
number of one of the hyperfine states, the resulting ground state can be imaged in a way 
that discriminates between them and the bubble structure may become visible. A high 
concentration of the denser particle species will accumulate at some point (s) in the trap. 
The division of the space between the BCS and the normal components is determined by 
the value of the gap thus, by studying its variation with the variation of the asymmetry, we 
can infer the existence of superfluidity and even the value of the gap, which is currently an 
outstanding problem. A better understanding of the surface tension of the interface between 
the two components is necessary to make this proposal fully quantitative. The qualitative 
arguments discussed above should be valid even if weak coupling BCS theory is not, as it is 
the case in the experiments with 6 Li [l7| . 
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